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A  RAY  THEORY  •‘OR  NONLINEAR  SHIP  WAVES 
AND  WAVE  RESISTANCE 


Hohyun  Ytm 

David  W.  Taylor  Naval  Ship  Research  and  Development  Center 
Be  the  sda  ,  Maryland  200h4,  U.S.A. 


Abstract 

Analytical  and  numerical  methods  for 
application  of  ray  theory  In  computing  ship 
waves  are  Investigated.  The  potentially 
important  role  of  ray  theory  in  analyses  of 
nonlinear  waves  and  wave  resistance  is  demon¬ 
strated.  The  reflection  of  ship  waves  from 
the  hull  boundary  Is  analyzed  here  for  the 
first  time. 

When  a  wave  crest  touches  the  ship  sur¬ 
face,  the  ray  exactly  follows  the  ship  surface. 
When  the  wave  crest  is  nearly  perpendicular  to 
the  ship  surface  thu  ray  is  reflected  many 
times  as  it  propagates  along  the  ship  surface. 
Many  rays  of  reflected  elementary  waves 
Intersect  each  other.  The  envelope  to  the  first 
reflected  rays  forms  a  line  lilte  a  shock 
front  which  borders  the  area  of  large  waves 
or  breaking  waves  near  the  ship. 

For  the  Wigloy  hull,  lay  paths,  wave 
phases,  and  directions  of  elementary  waves 
are  computed  by  the  ray  theory  and  a  method  of 
computing  wave  resistance  la  developed.  The 
wave  phase  is  compared  with  that  of  linear 
theory  as  a  function  uf  ship-beam  length 
ratio  to  identify  the  advancement  of  the  bow 
v  ave  phase  which  Influences  the  design  of  bow 
bulbs.  The  wave  resistance  of  the  Wigley  hull 
is  computed  using  the  amp'ltude  function  from 
Michell's  thin  ship  theory  and  compared  with 
values  of  Michell's  wave  resistance.  The  total 
wave  resistance  has  the  phase  of  hump  and  hollow 
shifted  considerably, 

lntroduct ion 


Significant  developments  in  the  ship  wave 
theory  have  been  made  in  recent  years.  These 
include  the  application  of  ray  theory'***  and 
the  experimental  discovery  of  a  phanomenon 
called  a  free  surface  shock  wave. 3 

Because  a  ship  is  not  thin  enough  to  apply 
both  the  thin  ship  theory  and  the  complicated 
free-surface  effect,  theoretical  development 
of  an  accurate  ship  wave  theory  has  been  slow. 
The  problem  should  be  evaluated  differently 
from  the  conventional  means.  Ray  theory  has 
been  used  in  geometrical  optica  or  for  waves 
having  small  wave  lengths.  Ureell^  used  the 


4A  complete  listing  of  references  is  given 
on  page  ! S 


theory  to  consider  wave  propagation  in  non- 
uniform  flow,  and  Shen,  Mayer,  and  Keller^  used 
It  to  investigate  water  waves  in  channels  and 
around  islands.  Recently  Keller*  developed  a 
ray  theory  for  ship  waves  and  pointed  out  that 
the  theory  could  supply  utaful  information 
about  the  waves  of  thick  ships  at  relatively 
.low  speeds.  However,  he  hed  difficulty  ob¬ 
taining  the  excitation  function  for  wave  ampli¬ 
tude  end  solved  only  e  thin-ehtp  problem.  Inul 
and  Kajitani'  applied  the  Uraell  ray  theory^  to 
ship  waves,  using  the  amplitude  function  from 
linear  theory. 


Because  the  ray  is  the  path  of  wava  energy, 
it  is  not  supposed  to  penetrate  the  ship  surface 
and  this  was  emphasized  by  Keller.  However, 
neither  Keller  nor  Inui  and  Kajltanl  consider 
nonpenetration  of  the  ray  seriously.  Recently 
Yim“  found  the  existence  of  rays  which  emanated 
from  the  ship  bow  and  reflected  from  the  ship 
surface, 


In  the  preaent  paper,  further  study  haB 
revealed  many  bow  rays  reflecting  from  the 
ship  surface,  cresting  an  ares  in  which  thase 
rsys  intersect  each  other.  The  envelope  to 
the  first  reflected  reve  forms  a  line  like  e 
shock  front  which  borders  the  area  of  large 
waves  or  breaking  waves  near  the  ship.  This 
will  be  referred  to  ss  the  second  caustic.  Thl6 
phenomenon  was  observed  by  Inul  et  el.*  in  the 
towing  tank  and  called  a  free  surface  shock 
wove  IFSSW) ,  ■  iu.  h  has  beta  done  concerning  the 

theoretical  in\ -  cigation  of  FSSW  by  introduc¬ 
ing  a  fictitious  depth  for  a  shallow  wafer  non- 
dlspersive  wave. 

The  ray  equation  and  the  ship  boundary 
condition  are  analyzed  further  to  show  the 
existence  of  reflecting  rays,  except  in  the 
case  of  a  flat  vedgelike  ship  surface.  The 
ray  path  is  vary  sensitive  to  the  initial  bound¬ 
ary  condition  near  the  bow  or  the  stern,  snt’ 
should  be  identified  at  downstream  infinity 
not  by  the  initial  condition.  The  conservation 
law  of  wava  energy  in  the  nonuniform  flow  is 
different  from  that  of  the  uniform  flow  due  to 
the  exchange  of  energy  with  the  local  flow. 8 
Therefore,  the  linear  wave  amplitude  aa  a 
function  of  initial  wava  angle  near  the  bow 
or  stern  is  meaningless  and  cannot  be  used  in 
the  ray  theory.  It  is  shown  that  the  amplitude 
function  for  the  ray  theory  matched  with  the 


t 


linear  theory  far  downstream  la  raaaonabla  and 
likely  to  produce  a  reasonable  result  as  In 
the  two  dimensional  theory.’ 

The  wave  phase  In  the  ray  theory  Is  ob¬ 
tained,  together  with  the  ray  oath,  Indepen¬ 
dent  trom  the  amplitude  function,  both  the 
phase  and  the  ray  path  are  quite  different  from 
the  prediction  of  linear  theory  near  the  ship. 
The  rays  far  downstream  are  straight  as  in  the 
linear  theory  yet  have  phases  different  from 
the  linear  theory;  they  are  parallel  to  the 
linear  ray  with  the  same  wave  angle  but  do 
not  coincide.  The  difference  in  wavs  phases 
is  the  main  factor  that  makes  the  wavs  resis¬ 
tance  different  from  the  linear  theory,  key 
paths  and  phase  difference  are  computed  for 
vurUnn.  parameters  of  the  Wlgley  Ship,  with  end 
without  a  bulbous  bow.  The  ray  paths  for 
different  drafts  and  different  beem-length 
ratio  of  the  Ulgley  hu’l  are  alightly  different, 
widening  the  wave  area  near  the  hull  for  the 
wider  beam  and/or  larger  draft.  However,  the 
phase  difference  Is  more  sensitive  to  the  beam- 
length  ratio  and/or  draft-length  ratio,  by 
always  advancing  the  linear  wave  phase. 

The  weve  resistance  of  the  Ulgley  hull  li 
computed  and  shown  to  have  e  considerable  shift 
of  phase  of  hump  and  hollow. 

The  most  Interesting  phenomenon  of  a  ahip 
with  a  bulbous  bow  ia  tha  reduction  of  slopes 
of  rays  and  the  second  caustic,  l.a.,  the 
larger  the  bulb,  the  greater  the  reduction. 

This  fact  was  observed  In  the  towing  tank.3 
There  exists  a  bulb  size  which  totally  elimi¬ 
nates  the  reflecting  rays.10  However,  the  phase 
dliference  due  to  the  bulb"  is  very  small 
showing  that  the  phase  difference,  which  has 
been  observed  in  the  towing  tank,  is  the  effect 
of  nonuniform  flow  caused  by  the  main  hull. 

Ray  Equations 


The  concept  of  ray  theory  In  ship  waves  la 
analogous  to  the  concepts  used  In  geometrical 
optics  and  in  geometric  acouatics.  A  ship  la 
considered  advancing  with  a  constant  velocity 
-li  which  lc  the  direction  of  the  negative  x 
axis  of  a  right  handed  rectangular  coordinate 
system  O-xyz  with  the  origin  0  at  tha  ahip  bow 
on  the  mean  free  surface,  t»0,  t  Is  positive 
upwards . 

First  the  phase  function  a(x,y,x)  it  de- 
f lned  so  that  the  equation  (a  •  constant)  ra- 
represents  the  wave  front  where  the  value  of  ■ 
is  the  optical  dlitance  fret  the  wave  source, 
c.g,,  the  ahip  bow.  Vhen  Keller*  developed 
his  ray  theory  of  shlpn  by  ex' ending  boundary 
conditions  and  a  solution,  which  should  aatlafy 
both  the  Laplace  equation  and  tha  boundrry 
conditions,  in  e  aeries  of  Froudo  numbs,  equaree 
F  ,  he  obtained: 


sz  -  -1  (S„  +  V$ 


(l) 


?•)*  at  a  •  0 


Sv  eliminating  S,  from  these  two  equations  he 
obtained  e  dispersion  relation 


<•,*  ♦  -,V 


(S,  +  * 
t  x 


4  s  ) 

y  y 


(2) 


where  4  is  the  double  model  potential.  When 
steady  state  motion  is  assumed  with  respect  to 
the  moving  coordinate  system,  0  -  xyr, 


When  the  angle  between  the  normal  n  to  the 
phase  curve  s  and  the 'axis  is  denoted  by  3, 

n  «  Icoaf)  a  jslne  (3) 

Then  the  wave  number  vector  Is  defined  by 


k  =  k,i  +  k2j  ;  sxl  +  syj 


•  nk  »  Ik  cob9  +  jk  sine  it  :  *  0 
From  equations  (2)  and  (4) 

k.r _ ! _ i2 

fu  coafl  +  v  sine] 


(4) 


(5) 


where 


-4  ■  u  and  - p  ■  v 

x  y 

These  results  are  an  approximation  within 
the  order  of  F*,  end  the  phase  function  and  its 
ralated  equations  are  all  limited  to  their 
values  at  2  >0.  Thus,  from  now  on,  unless 
otherwise  mentioned,  all  the  physical  values 
are  at  t  »  0.  The  ray  equation  of  ship  waves 
la  obtained  from  the  irrotat tonality  of  the 
wave  number  vector 


■Ik, 


3k 

— i  •  0 
3y 


Frota  Equal  long  (A)  through  (6), 

{2  »in0  (u  j  im0  -  v  cost)) 

4  coed  (u  cose  +  v  sine)} 

4  {2  cosb  (— u  tin?  4  v  cose) 

4  ainb  (u  cose  4  v  sine)}  7^ 

'  'y 

*  2  line  (cose  +  gtnp  il; 

-  2  cos6  (cos6  ~  +  sln9  7—) 

3y  V 


(6) 


(?) 
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This  la  the  ray  i-quatlon  which  can  b« 
solved  by  the  method  of  characteristics,  and 
Is  equivalent  to  simultaneous  ot dinary  differ¬ 
ential  equations. 

dy  (v-  S  ainO  (u  cost*  +  v  sinO)}  ^ 

dx  {u-  S  co»n  (u  cosu  +  v  atari)} 

i)  t1  3u  \v 

“(2  slno  (cost'  *jj  +  slnB  — ■) 

-  2  COP0  (coaO  ff  +  nine  |£>}/<2  sinB  (u  slnB (9) 

-  v  co*8)  e  cosO  (u  cosO  +  v  sint))l 

Here  u  cosB  s  v  sin  O  Is  the  velocity  component 
normal  to  the  phase  curve  of  the  flov  relative 
to  the  ship.  Because  the  wave  la  stationary 
relative  to  the  ship,  the  phase  velocity 
through  the  water  surface  should  be 


through  the  stagnation  point  where  u  end  v  are 
the  velocity  components  of  the  total  velocity 
caused  by  the  double  model  source  distribution 
m  and  the  uniform  flow  relative  to  the  ship. 

In  the  linear  ship  wave*  theory  a  smooth 
source  distinction  produces  two  systems  of 
regular  ship  waves:  the  bow  and  the  stern 
waves  starting  Torn  the  bow  and  stern, 
respectively,  represented"  as 

v/e 

4(:;,y)  »  J  A.(t3)  exp  {l  Sj  (x,y,0)}  dO  (14) 
-n/2 

where 

Sj  »  kQ  sec^O  {(xmXp)  co»8  +  y  sine}  (15) 
*  the  x  coordinate  of  the  bow  or  ster-* 


-  u  cosB  -  v  slnB 


The  group  velocity  la  one-half  of  tha  phase 
velocity.  Thus  the  ray  direction  in  Equation 
(B)  is  along  the  resultant  of  the  group  velo¬ 
city  taken  normal  to  the  phase  curve  snd  of  the 
velocity  of  the  basic  flow  as  Ursell  has  shown, ^ 
Because  the  wave  energy  is  propagated  at  the 
group  velocity,  the  ray  path  la  lntcrprettd  as 
the  path  of  energy  relative  to  the  skip.  This 
can  be  obtained  by  solving  Equations  (6)  and 
(9)  with  tha  proper  Jritial  condition.  The 
phase  a  can  be  obtained  from  aquations  (4)  by 

a  -Jk  dr 


a r.  In  potential  thaory;  a(x,y)  la  a  function  of 
(x.y)  but  la  unrelated  to  the  Integration  oath. 
However , 


Is  -  k  cost*  dx  +  k  sin 6  dy  (10) 


can  he  solve,  .ogether  with  tha  ray  equations 
along  tha  ray  path. 


Raya  of  Ship  Waves  and  Linear  Thaory 

To  invest igata  tha  path  of  a  ray  of  a 
ship  wave,  a  ahlp,  rapraaantad  by  a  double 
modal  source  distribution  m  (x.y)  on 
y  •  o,  h^x^-h,  is  coc.aldared,  Although  tha 
linear  relation  between  tha  source  strength 
m  and  the  ship  surface 


y  -  If (x,») 


Is 


no 


(12) 


the  actual  double  model  ahlp  body  straamilna 
should  be  obtained  by  solving 


dy  ^  v 

3x  u 


(13) 


g  »  acceleiation  due  to  gravity 

/.(B)  ■  amplitude  function  which  la  a 

function  of  source  distribution  m 

The  regular  wave  4  Is  the  solution  of 
llnaar  ahlp  wave  thaory  far  from  tha  chip." 
Actually,  It  la  aaay  to  see. that  tha  exponential 
function  satisfies  loth  tha  laplaca  aquation 
and  tha  llnaar  free-surfAce  boundary  condition 
for  any  value  of  xi ,  Havelock  Interpreted 
Equation  (14)  by  discussing  tha  Integrand  aa 
elementary  waves, 12,13  l.a.,  tha  regular  waves 
are  aggregates  of  elementary  waves  starting  from 
tha  bow  and  atam  of  a  ship.  Tha  normal 
direction  of  each  elementary  wava  craat  la 
n  “  X  odsB  +  J  alnB. 

Becauaa  the  local  disturbance  of  the  double 
model  decays  rapidly  away  from  the  ship,  even 
in  nonlinear  theory  in  the  fer  field,  the 
regulir  wave  should  b.  of  the  form  of  Equation 
(14)  with  possibly  different  values  of  xi  and 
A(B) .  When  tha  integral  of  4  is  evaluated  by 
tha  method  of  tha  stationary  phaaa1*  by  taking 
roots  of 


|f  (x.y ,o)  -  o 


(16) 


2  tan^B  +  f  tanB  +  1  •  0  (17) 

in  general  two  values  of  B  are  obtained  for  a 
given  value  of  each  x  and  y,  aatlsfying, 

>  B*1  W8) 
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itiul  Hu*  vaiue  o t  the  Integral  of  Equation  (14) 
lor  any  x.y  In  x/|y|2*8*3  can  be  evaluated  by 
llu'  stat  l  unary  i>hase  method  as  the  sum  of  two 
waves:  transverse  waves  (0  de^  <  { 6 (  <  .15  deg) 
and  divergent  waves  (35  deg  <  T 0 }  <  90  deg). 

Except  for  the  amplitude  function  exactly 
ilu*  same  result  for  the  relation  of  x.y  and  0 
.is  In  the  linear  theory  can  be  obtained  from 
Equation  (8)  by  substituting  u  *  l  and  v  M  0, 
Thus,  this  means  chat  the  energy  of  each 
elementary  wave  propagates  on  the  uniform  flow 
along  a  straight  line  ray.  The  rays  are  con- 
t  ined  u\  x/jy|  i  8s  both  of  the  elementary 
waves  at  u  ■  0  and  0-90  deg  correspond  to  the 
vuy  %•  /\  *•  0,  and  fl  ■  35  deg  corresponds  with 
tin*  tav  x/|yj  •  8».  One  ray  between  these  two 
lays  corresponds  to  two  elementary  waves  where 
one  is  trausversel  and  the  other  is  divergent. 

When  the  velocity  is  not  uniform  due  to  the 
t low  perturbation  caused  by  a  ship,  the  ray  is 
not  straight  but  curved  near  the  ship  as 
shown  in  Figure  1.  it  is  known  that  wave  energy 
t  lux  divided  by  the  relative  frequency  with 
lospect  to  the  coordinates  for  which  the  fluid 
vc ! sv  it y  is  rero.  is  conserved  along  the 
cuived  iuy  lube.  If  the  coordinates  are  fixed 
in  space,  then  the  ship  waves  are  unsteady 
illative  to  thi*  fixed  coordinates  and  the 
;  icqucm  y  is 


s ^  -  U  k  cosu 


(19) 


Figure  l  -  Ray  Pa:hs  for  a  Wigley  Hull 
(b  -  0.1.  h  -  0.0825) 


because  only  the  transformation  ol  coordinates 
(x.y)  to  (x  -  Ut i  y)  Is  needed  to  get  the  un¬ 
steady  flow.  Strictly  speaking  the  relative 
frequency  Is 

■  k  u  co»0  +  k  v  sinu  -  k  U  cos" 

but  tor  our  discussion,  the  approx  Invite  value 
serves  the  purpose.  This  means  that  for  a 
uniform  flow  the  wave  energy  is  constant 
along  the  ray  because  for  a  unltorm  flow  u  - 
V  -  0.  K  -  g/(l)2  cosAi)  ,  and  11  Is  constant 
along  the  straight  tay.  However,  the  wave 
energy  la  not  constant  along  the  curved  ray  In 
nonunifont  flow,  but  la  dependent  upon  the 
local  velocity  components  and  0  because  the 
relative  frequency  la  a  function  of  the  local 
velocity  and  9  as  In  Equations  (5)  snd  (19)  and 
u  changes  along  the  curved  ray.  The  wave 
energy  will  change  considerably  alonR  the 
curved  ray  near  the  bow  or  stern  because  u  and 
v  change  to  lero  at  the  wave  source,  the 
stagnation  point,  while  In  the  far  field  lv 
will  be  constant  along  the  ray  as  In  linear 
theory.  Near  the  stagnation  point,  the  wave 
number  Increases  according  to  Equation  (5)  and 
the  wave  energy  also  Increases  due  to  the  energy 
conservation  law.®  Thus,  the  wave  may  break  near 
the  stagnation  point. 

Becauae  wav*  amplitude  Is  so  difficult  to 
obtain  from  the  ray  theory^,  the  linear  theory 
haa  been  conaidered  as  an  approximation.!  A 
good  reault  was  obtained  for  a  two  dimensional 
flow  problem. ^  However,  extreme  care  Is  needed 
In  the  three  dimensional  theory  because  0  chang¬ 
es  by  a  large  amount  along  the  ray,  neat  the 
stagnation  point,  and  the  wave  energy  depends 
upon  0.  The  matching  amplitude  function  of  ray 
theory  with  the  values  from  linear  theory  ahuulJ 
be  done  at  the  far  field  where  6  end  the  direct¬ 
ion  of  the  ray  ate,  respectively ,  identical  for 
both  cases  for  each  elementary  wave.  In  addi¬ 
tion,  the  Initial  condition  has  to  be  taken  in 
the  neighborhood,  but  not  exactly  at,  the  stag- 
national  point  becauae,  although  the  stagnation 
point  is  the  weve  source,  at  the  stagnation 
point  the  ray  equations  are  Indeterminant ,  a 
and  v  being  tero.  However,  because  of  the  large 
change  of  u  and  v  near  the  stagnation  point,  the 
ray  path  1b  very  sensitive  to  the  Initial  values 
of  x,  y  and  9;  the  nearer  to  the  stagnation 
point,  the  more  sensitive.  Therefore,  the  iden¬ 
tification  of  ray  paths  should  :e  correlated 
with  the  values  of  s  at  inf  Ini' y.  Then  all  the 
ray  paths  can  be  properly  and  uniquely  identified 
by 
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Since  perturbations  of  the  ship  decay  . 

rapidly  away  from  the  ship,  0  also  rapidly  j 

approaches  9».  The  relation  between  the  initial  1 

value  end  the  value  at  Infinity  of  9  has  little  I 

meaning,  although  it  was  misunderstood  before"6 
because  9  changes  very  rapidly  ne*\  the  stag¬ 
nation  point.  J 
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Rays  of  Ship  Waves  and  Ship  Boundaries 

Ship  waves  created  by  a  smooth  ship  hull 
propagate  as  regular  bow  and  stern  waves  from 
the  bow  and  stern  stagnation  points  to  infinity 
along  rays.  Because  a  ray  carries  the  wave 
energy  it  cannot  penetrate  the  ship  surface.  If 
the  linear  tree-surface  condition  is  considered 
with  the  exact  hull  boundary  condition!  as  has 
been  popular  in  recent  ship  wave  analysis^, 
the  straight  rays  which  pass  through  the  hull 
boundary  have  to  be  considered.  For  a  ray 
theory  with  the  exact  hull  boundary  condition, 
the  ray  is  not  allowed  to  penetrate  the  ahip. 

In  fact,  when  the  initial  wave  crest  touches 
the  ship  boundary  it  can  be  proved  that  the 
ray  of  such  a  vave  graces  along  the  ship 
boundary  without  penetrating  the  ship  boundary. 


When  the  wave  crest  is  perpendicular  to 
the  ship  hull 

^  *  tanO  (26 ) 

If  Equation  (2b)  is  inserted  into  Equation  (8) 

&  *  t,nu  *  z  (27> 

The  ray  also  touches  the  ship  Initially. 
However,  Equation  (9)  Is  not  compatible  with 
Equation  (26)  on  the  hull.  This  can  be  proved 
In  a  similar  way  as  follows: 

Differentiating  Equation  (26)  with  respect 
to  x  along  the  ship  hull 


In  Equation  (8)  when  the  wave  creat 
touches  the  hull 

u  cost?  +  v  sine  »  0  (20) 

because  t?  is  the  angle  between  the  normal  to 
the  crest  and  to  the^axis,  and  the  velocity 
normal  to  the  ship  hull  is  zero  on  the  hull. 
Equation  (8)  then  becomes 


dy  ^  v 
dx  u 


de  _  * xx 

dx  PTI 


(28) 


Inserting  Equation  (26)  Into  Equation  (9)  gives 


***  u  (1  +  —  tan3) 

u 

From  Equations  (25) i  (27),  and  (29),  noting 


showing  that  the  ray  touches  the  ship  hull 
streamline  from  Equation  (13). 


3u  jv 
Is  "  3x 


When  the  wave  crest  touches  the  hull, 
from  Equations  (11)  and  (13) 


v 

u 


(21) 


and  from  Equations  (20)  and  (21) 

tanO  -  -p-  (22) 

X 

Differentiating  Equation  (22)  with  respect  to 
x  along  the  ship  hull 


dx 


(23) 


x 

However,  inserting  Equation  (20)  into  Equation 
(9)  yields 


do 

n 


..  3u  ^  iv  2.  3u  3v 

cott)  —  +  —  -  cot  6  r-  -  cot9  t — 

3x  3x _ 3y _ 3y  (24) 

u(l  -  —  cot9) 


—  -  -o  ■---  (30) 

dx  f  +  1 
x 

Equations  (28)  and  (30)  are  compatible  only 
when  £xx  ■  O,  or  fx  •  constant , 

That  is,  only  when  the  ship  is  a  flat  plate 
does  the  ray  of  the  vave.  whose  crest  is  perpen¬ 
dicular  to  the  ship,  follow  the  ship  boundary. 
This  means  that  when  the  ship  bow  is  like  a 
wedge,  the  ray  of  the  bow  wave,  whose  crest  Is 
perpendicular  to  the  wedge  surface,  initially 
follows  the  wedge  surface, 

When  the  ray  equations  are  solved  numer¬ 
ically  by  the  Runge-Kutta  method,  with  initial 
values  near  the  bow  stagnation  point,  the  ray 
touches  and  follows  the  ship  boundary  at 

81  -  -  3  +  C.  (31) 

where  a  is  the  half  antrance  angle  of  the  ship 
bow  and  8j  denotes  tha  Initial  value  of  6. 


Differentiating  Equation  (21)  with  reipec*  to 
x  along  the  hull  yields 

v  iH  +  2^  v  3u  3y  3v  3y 

3  v  _  ~  u  Jx  3x  ~  u  3y  3x  3y  3x  (25) 

OX  U  SJ 

From  Equations  (20),  (21),  (22),  and  (25) 
it  can  be  shown  that  Equations  (23)  and  (24) 
are  equivalent.  That  is,  when  the  wave  crest 
touches  th>-  ship  boundary,  the  ray  equations 
and  the  ship  vull  streamline  equations  are 
equivalent . 


When  Si  increases  from  this  value  the  ray 
aovas  gradually  away  from  the  ship  as  shows  in 
Figures  1  through  4.  The  ray*  are  curved  near 
the  ship  but  at  far  downstream  they  are  straight 
and  the  ray  direction  becomes  exactly  the  same 
function  of  8  as  the  linear  theory.  Thus  at 
infinity  the  rays  are  inside  |dx/dy|  *  o1. 
However,  when  S»  approaches  sero,  the  curved 
ray  near  the  ship  gradually  approaches  the 
ship,  and  eventually  crosses  the  ahip  boundary, 
as  ir  Figures  1  through  4.  Here  the  wave 
reflection  should  be  considered  to  prevent  the 
ray  penetration  of  the  ship  hull. 
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tgure  2  -  Rav  Paths  for  a  Wlgley  Huli 
(b  -  0.2,  h  -  0.0625) 


tit 


X 

Figure  3  -  Ray  Paths  for  a  Wlgley  Hull 
(b  -  0.2,  h  -  0.03) 


Kay  Reflection 

There  are  no  waves  coming  from  places 
other  than  the  ship  bow  or  the  ship  stern.  In 
the  ray  theory,  the  flow  field  perturbed  by  the 
ship  deflects  the  ray  path  starting  from  the 
ship  bow  Loward  the  ship.  Thus,  some  rays 
of  bow  waves  impinge  on  the  ship  hull. 

Because  the  ray  theory  is  for  small  Froude 
timbers  pp.d  th»  wave  phenomenon  is  considered 
only  on  the  free  su'face,  it  is  reasonable  to 
consider  only  reflected  waves,  as  in  geometrical 
optics,  neglecting  transmitted  waves  when  the 
oncoming  waves  impinge  on  the  ship  hull. 

Let  *  e  ray  at  the  wave  angle  0  on  the 
ship  bounuary  (x,y)  be  reflected  to  er  and 
fx  -  tan0o  as  in  Figure  5. 

Then , 


or 

e  -  26  -  0  (32) 

r  o 


Figure  5  -  Angles  of  Waves  Reflecting  on 
the  Ship  Hull 

Whenever  the  ray  intersects  the  ship 
boundary  at  (x,y)  with  angle  u,  then  the  value 
of  O  at  (x,y)  will  be  changed  to  t)r  "  20o  -  0, 
and  x,y*0  are  continuously  calculated  by  the 
Runge-Kutta's  method.  Then  the  ray  will  reflect 
as  in  Figures  l  through  4.  Here  if  Q0  is  zero, 
it  is  easily  seen  in  Equation  (8)  that  the 
impinging  ray  angle  (J)  changes  to  -(H)  for  the 
reflected  ray  angle.  This  fact  can  be  easily 
shown  to  be  true  even  for  0y  *  0  by  just  the 
rotation  f  the  coordinate  system. 


Nunver  cal  Experiments  of  Ray  Paths  and 
Free  Surface  Shock  Waves  ~~ 

Becaus '  the  ray  equations  an  be  solved 
only  numeric,  lly,  careful  numerical  experiments 
with  the  rey  equations  may  give  valuable  in¬ 
formation.  For  simplicity,  the  Wigley  hull 
source  distribution  with 


2tt  dx 


where 

f j  -  2irb  (1  -  x2)  {l  -  (|)2)  (33) 

is  considered  for  various  numbers  of  b  and  h 
which  are  related  to  the  hull  beam  and  draft, 
respectively.  The  sctu  1  hull  shape  corre¬ 


sponding  to  the  source  Equation  Oil)*  is  ob¬ 
tained  by  plotting  the  body  streamline  passing 
through  stagnation  points  as  the  solution  of 

dy  v 

dx  u 

where  u  and  v  corresponding  to  Equation  (32) , 
are  shown  in  Appendix  A  together  with  u^,  u  , 
and  vy  for  the  ray  equation.  Equations  (8)\ 

(9),  and  (lOj,  together  with  the  screamline 
equation,  are  solved  by  the  Runge-Kutta  method 
with  initial  conditions  (x,y,0)  near  the  stag¬ 
nation  points  with  various  values  of  0. 

Many  ray  paths  both  reflecting  and  non¬ 
reflecting  from  the  surfaces  of  various  Wigley 
hulls  are  shown  in  Figures  1  through  4,  These 
paths  were  computed  by  a  high  speed  Burroughs 
computer  at  David  Taylor  Naval  Ship  Research 
and  Development  Center.  The  reflection  condi¬ 
tion  is  incorporated  in  the  high-speed  computa¬ 
tion  with  a  routine  to  find  the  intersection  of 
the  tay  and  the  ship  boundary  which  is  pre¬ 
calculated  and  saved  in  the  memory.  The  step 
sizes  of  integrations  and  interpolation  were 
determined  after  many  numerical  tests,  and  the 
sh'  vn  results  are  considered  to  be  reasonably 
accurate.  For  a  given  initial  condition,  the 
solution  is  stable  and  converges  well. 

In  Figures  1  through  4,  some  common  fea¬ 
tures  of  rays  can  be  drawn  as  follows.  The 
rays  in  -  z/2  <  8®  <  0  far  behind  the  ship 
oehave  like  rays  of  linear  theory  except  wave 
phases  are  advanced  in  the  ray  theory  and  those 
near  0®  ■  0  are  rays  propagating  from  the  ship 
bow  and  reflecting  *roro  those  of  ship  hull. 

The  rays  near  the  ship  are  very  different  from 
the  linear  theory  as  Inui  and  Kajitani1  pointed 
out.  The  curved  rays  from  the  bow  have  a  far 
larger  slope  than  those  of  linear  theory  and 
the  phase  of  each  ray  is  considerably  advanced. 
The  magnitude  of  the  phase  difference  is  more 
sensitive  to  the  beam-length  ratio  and  the 
draft-length  ratio  than  the  magnitude  of  the 
ray  slope. 

When  6*,  «  35  deg,  the  ray  will  be  the 
outermost  ray  and  the  ray  angle  at  ®  will  be 
approximately  tan"1  8'^  as  in  the  linear  theory 
and  there  is  the  corresponding  initial  value  of 
0  or  Bi  near  the  bow,  or  the  origin.  However, 
the  Oi  which  is  corresponding  to  a  single  value 
of  0®  is  very  sensitive  to  small  changes  of  x 
and  y  near  the  origin.  At  a  fixed  point  near 
the  origin  there  exists  a  unique  correspondence 
between  6^  and  6oo. 

When  from  the  @i  which  is  corresponding  to 
8*  «  35  deg,  the  initial  value  of  0^  increases, 
6«>  also  increases  and  the  ray  angle  decreases. 

In  general,  when  Oi  •  0,  0®  is  still  a  negative 
value.  When  0i  increases  further,  approaches 
zero  and  the  ray  pass  is  very  close  to  the 
stern.  At  this  point,  there  exists  a  certain 
increment  of  0^  which  makes  the  ray  barely 
touch  the  ship  stern,  at  3^  *  0io.  When  Oi 
slightly  increases  from  0io,  the  ray  reflects 
from  the  ship  hull  ne«r  the  stern.  With  the 
incrament  of  ©i  the  reflection  point  moves 
toward  the  bow.  When  0  -  8^  the  ray  once 
reflected  touches  the  stern  again.  When  0i 
increases  further  from  0il,  the  ray  reflects 
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uwi'  iron*  the  ship  hull.  In  lhi:>  way,  further 
uni'MiH-nt  ol  si  j  makes  the  rav  reflect  Iron)  the 
^.hip  1ml  I  three,  four  ...  times.  However,  at 
!  roar  the  value  of  i,  the  ray  trios  to  pene¬ 
trate  the  ship  hull  at  the  starting  point  of 
the  bow.  This  cannot  be  allowed  because  this 
kind  of  ray  should  come  from  outside  of  the 
slu |  .  l.ol  the  border  point  of  Oj  be  0^u. 

Ibis  means  that  rays  of  initial  value  of  0 
between  ■■ll,  <  t! i  <  u iu  reflect  from  the  ship 
hull.  Ai.  i:>  clear  in  Figures  1  through  4, 
m  genera),  all  the  rays  before  reflection  do 
not  intersect  each  other,  however,  reflected 
ravs  Lntersect  other  rays.  The  once  reflected 
rays  intersect  not  only  with  each  ray  once 
reflected  from  ship  boundary  potnts  close  to 
each  other,  but  also  with  af  least  one  ray 
be 1  ore  re f 1 ec t ion . 

In  the  stationary  phase,  each  ray  has  an 
amplitude.  Likewise  in  the  ray  theory,  each 
lay  carries  it's  energy.  The  reflected  ray  may 
have  approximately  the  sane  energy  as  the  ray 
at  vj  *  v io  or  0*  H  o  where  the  amplitude 
function  of  the  linear  theory  is  in  general 
more  significant  than  amplitudes  of  the  other 
values  of  Uo..  Because  the  phase  must  be 
approximately  close  to  each  other  for  the  waves 
near  •'.»  -  0,  the  wave  height  of  the  once  re- 
I  levied  lay  may  be  close  to  three  Limes  that 
ol  the  transversal  wave  for  'Jo,  *  0.  When  the 
envelope  of  the  once  reflected  waves  is  drawn, 
i he  domain  bounded  by  the  ship  surface  and 
the  envelope,  denoted  by  Dm,  must  be  distinctly 
d 1 1 1 e rent  from  other  domains  because  in  Dm 
there  are  not  only  once  reflected  rays  hut  also 
twic  or  multi-reflected  rays  on  which  more  than 
three  tei levied  rays  intersect  by  an  argument 
similar  to  that  used  for  the  once  reflected 
rays.  The  envelope  of  the  once  reflected  rays 
behaves  like  the  shock  front  which  was  observed 


In  general,  a  line  is  called  a  caustic  when 
-n  side  of  the  line  one  can  find  a  continuous 
•  li  a  t  ibut  ion  .*t  rays  but  not  on  the  other  side, 
and  along  the  caustic  the  wave  slope  is  found 
io  be  large.  l‘he  wave  near  the  rav  angle  t‘r»-  « 
deg  is  the  caustic,  and  the  wave  height  near 
the  caustic  far  downstream  can  be  obtained  by 
an  application  of  the  Airy  function  in  the 
linear  theory.  The  envelope  of  the  once  re- 
tinted  lavs  tn.iv  be  a  kind  of  caustic  formed 
bv  the  let  tact ed  bow  wave  rays  due  to  the  non- 
nmtoim  flow  perturbed  by  the  ship.  Shen , 

M» vei ,  and  KellerJ  studied  such  caustics  caused 
bv  the  sloping  beach  of  channels  and  around 
islands.  Thus,  the  additional  caustic  of  ship 
waves  nuiv  be  called  the  second  caustic  of  ship 
waves,  ami  should  not  be  confused  with  the  first 
vauHt ir  which  is  the  known  caustic  at  0^  *  35 
•leg. 

:.n  oiul  t  aust  i c  of  Ship  Waves  of  Various  Ships 


For  the  Wigley  hull,  several  different 
values  of  parameters  b  and  h  were  taken  to 
find  their  effect  on  the  second  caustic.  In 
addition,  the  effect  of  a  bulbous  bow  on  the 
second  caustic  was  considered.  The  most 
.1 1st  iugulstiable  physical  characteristic  of  the 
se«  omJ  oaiist  iv  is  its  distance  from  the  ship 
hull.  This  is  related  to  the  distance  of 


reflected  rays  from  the  ship  hull.  If  the 
number  of  reflection  rays  Increases,  or  If 
increases  from  Bi0,  the  maximum  distance  be¬ 
tween  the  ray  and  the  ship  hull  decreases. 

The  distance,  before  or  after  the  ray  reflection, 
approximately  behaves  like  a  sine  curve.  The 
maximur  distance  between  the  ray  before  the 
first  reflection  and  the  ship  hull  divided 
by  the  x  coordinate  of  the  point  of  the  first 
reflection  a/xj  is  plotted  in  Figure  6  for 
various  ships.  The  value  of  a/x[  for  different 
values  of  0i  are  approximately  the  same  for  a 
given  hull  and  are  related  to  the  area  between 
Lhe  second  caustic  and  tha  ship  hull  whe.re  there 
may  be  breaking  waves  or  turbulent  waves.  Thus, 
if  the  area  ia  large,  viscous  dissipation  of 
energy  becomes  large.  Accordingly,  the  measured 
momentum  loss  behind  the  ship  for  the  breaking 
waves  becomes  large , 

The  values  of  a/xj,  Increase  with  increas¬ 
ing  beam-length  ratio.  However,  the  most 
interesting  part  is  the  effect  of  the  bulbous 
bow. 10  When  the  bulb  size  is  Increased  or  the 
doublet  strength  is  increased  the  curvature  of 
the  ray  near  the  bow  becomes  leas,  although  the 
streamline  near  the  bow  is  such  that  the  en¬ 
trance  angle  is  slightly  large  The  values  of 
a/x]  decrease  with  increasing  bulb  size,  and 
eventually  the  ray  for  6W  K  0  cannot  propagate 
without  penetrating  the  ship  hull  at  the  begin¬ 
ning.  That  is,  there  is  no  reflecting  ray 
coming  out  of  the  bow  with  a  bulb  of  the  proper 
size,  as  shown  in  Figure  7. 


Figure  6  -  Width  of  the  Second  Caustic 
a/xj  for  Wlgley  Hulls 
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Figure  7  -  Rays  of  the  Wigley  Ship  (b  »  0.2,  h  »  0.0625)  with  Bulb 
(radius,,  r^  •  Q.02l4h>  depth*  Zj  ■  0.5h) 


These  phenomena  associated  with  the  second 
caustic  are  exactly  the  same  as  those  of  the 
"free  surface  shock  wave"  which  was  observed 
in  the  towing  tank. 

Finally  it  should  be  pointed  out  that  the 
ray  of  the  wave  whose  crest  is  perpendicular 
to  the  hull  surface  follows  exactly  along  the 
flat  surface  as  was  proven  by  Equations  (26) 
and  (30)  and  the  following  paragraphs.  This 
means  that  rays  originating  from  the  vertex 
of  a  wedge  vary  likely  nevr r  reflect  on  the 
wedge  surface.  This  is  because  toys  near  the 
ray  which  follows  the  ship  surface  did  not 
intersect  each  other  near  t*i  -  -*/2  4-  u,  or 
m  -  "/2.  Rays  near  0W  X  0  impinge  on  the 
ship  surface  due  to  the  effect  of  the  water-* 
line  curvature  of  the  ship.  Therefore,  if  a 
ship  has  a  wedge  bow,  the  second  caustic  must 
be  found  near  or  behind  the  shoulder.  Be¬ 
cause  the  first  caustic  near  the  bow  may  be 
very  prominent  and  both  waves  near  the  first 
caustic  and  waves  near  the  stagnation  point 
break,  careful  experimental  analysis  and  more 
theoretical  study  of  the  bow  near  field  may  be 
needed . 


Wave  Amplitude  and  Phase 


A(0)  -  P(0)  +  i  QUO  (34) 


may  be  taken  as  an  approximation  from  the  linear 
theory  but  the  phase  difference  xi  must  be  ob¬ 
tained  from  matching  with  the  near  field.  If 
the  near  field  is  also  expressed  by  the  linear 
theory 

o  1 


P  +  i  Q  »  n  kQ  /  / 


iiidz  f  e 


-h  o 

kQZ  sec  0  x  aecft 

sec3e  e  0 


(35) 


where  fx  is  the  derivative  of  Equation  (11)  with 
respect  to  x.  When  the  inner  integrand  of 
Equation  (35)  is  integrated  with  respect  to  x, 
the  value  with  the  limit  x  ■  l  will  form  stem 
waves  and  the  value  with  the  limit  x  ■  0  will 
form  bow  waves. H  Then  the  bow  waves  can  be 
represented  by 

tt/2 

■  f  \  cxp  8i  (W 

-n/2 


Because  the  perturbation  due  to  a  ship, 
or  both  regular  waves  snd  local  disturbances, 
decays  at  far  field,  the  linear  theory  oust 
hold  in  the  far  field.  In  particular,  the 
wave  resistance  can  be  calculated  from  the 
energy  passing  through  a  vertical  plane  xa 
constant  far  downstream;  the  linear  theory 
which  is  prcperly  matched  to  the  near  field  of 
the  ship  will  be  used  for  calculating  the  wave 
resistance.  As  explained  in  Equation  (14)  and 
the  following  paragraphs,  the  expression  for 
regular  waves  far  downstream  is  known  to  be  of 
the  form  of  Equation  (14)  where  the  amplitude 
function 


where 

s  <e>  -  pb  +  1  Qb  (37) 

When  the  phase  s  is  computed  {rum  Equation 
(10)  along  with  the  ray  path  from  Equations 
(8)  and  (9)  considering  that  s  »  0  at  the  bow 
near  the  origin,  there  are  two  results  different 
from  those  of  linear  theory:  (1)  the  ray  path 
is  deflected  as  if  the  elementary  wave  of  the 
linear  theory  started  from  (xj,0)  not  from  the 
origin,  end  (2)  the  phase  change  denoted  by 
As  should  be  considered,  that  is,  the 
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equivalent  linear  e  leimuiUiry  wave  m-~.y  be 
written  as 


Wave  Resistance 


.  2 

A,  ('•)  exp  [i  k  .sec  f> 

°  (38) 

|(x-Xj)  cost)  +  y  slno)  +  i  As] 

where  xj  Is  obtained  as  an  intersection  of  the 
tangent  to  the  ray  at  •»*  and  the  x  axis  and 

k^.'.s  *  V^s  -  k}  sec“M  |(x-Xj)  cos'1  +  y  sinoj 

(39) 

The  value  of 


if  all  the  elementary  waves  are  assumed  to 
be  propagated  without  reflection,  the  amplitude 
function  and  the  phase  difference  studied  in 
the  previous  sections  will  supply  enough  infor¬ 
mation  for  the  calculation  of  wave  resistance. 
Because  at  far  downstream,  the  wave  height 
may  be  considered  linear,  the  Havelock  wave 
resistance  formula^  may  be  used  for  the  waves 
represented  by  Equation?  (36)  thr >ugh  (AO), 
considering  that  the  wave  with  the  changed 
phase  82^  (0)  has  the  amplitude 

(P.  +  1  Q  )  eik°e2b  (b) 

□  b 


s  ,  (u)  *-  As  -  yj  seed  (AO) 

can  be  obtained  at  any  point  along  the  ray. 

In  general,  xj  is  negative  and  s2  ('0  is 
positive  meaning  that  the  bow  wave  phase  in  the 
ray  theory  is  larger  or  more  advanced  than  the 
phase  of  the  linear  theory.  This  fact  has  long 
been  observed  in  experiments  in  towing  tanks. 
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e/2  ITL 


(Al) 


ik  s9.  (0)  , 

|  (P^  +  iQ^)  e  (  cosJudh 


The  advancement  of  wave  phase  is  computed 
for  various  ships  and  the  values  of  s 2  and  xj 
at  x  -  2  are  shown  in  Figures  8  through  10. 

When  the  beam- length  ratio  Increases  and/or 
the  draft-length  ratio  increases,  the  values 
of  s;>  Increase  and  the  values  of  Xj  decrease 
for  all  valuer  cf  0U .  As  compared  with  the 
increment  of  the  slopes  of  rays  near  the  ship, 
the  increment  of  the  phase  angle  is  more 
sensitive  to  the  beam-  and/or  draf t-to-length 
ratios. 

The  most  interesting  phenomenon  about  the 
phase  difference  is  in  regard  to  the  bulbous 
bow. 10  That  is,  the  phase  differences  for  hulls 
with  and  without  bulbs  are  almost  the  same  even 
with  a  considerably  larger  bulb.  In  the  past, 
because  of  the  observed  phase  difference  of 
ship  waves,  the  bulb  was  located  far  forward 
to  obtain  good  bow  wave  cancellation. ^ 

According  to  the  present  analysis,  if  there  is 
no  other  reason,  the  bulb  position  need  not  be 
far  iorvatd.  Because  the  nonuniform  flow  cre¬ 
ated  by  t He  ship  is  much  more  significant  than 
that,  of  the  bulb,  as  far  as  phase  change  is 
concerned,  both  the  ship  bow  waves  (in  general, 
positive  sine  waves)  and  the  bulb  waves 
(negative  sine  waves)  propagate  through  the 
sHrae  region  and  cancel  each  other. 

As  for  the  amplitude  function,  although 
it  was  shown  by  Doctors  and  Dagan^  that  iny 
theory  produced  the  best  result  for  a  two- 
dimensional  submerged  body  even  though  they 
used  a  linear  amplitude  function,  the  surface 
piercing  three-dimensional  case  may  be  quite 
different.  The  amplitude  function  is  mainly 
related  to  the  singularity  strength  which  sat¬ 
isfies  the  ship  hull  boundary  condition  and 
some  improvement  might  result  by  considering  the 
sheltering  effect.  However,  in  the  present 
study,  the  linear  amplitude  function  is  used  to 
simplify  the  problem,  showing  the  effect  of 
curved  rays. 


>r  in  Srettensky’s  formula 
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Figure  10  -  S,a»J.  and  ~x  of  a  Wigley  Hull  (b  *  0.2.  h  »  0.0625)  with  and  without 
Bulb  (tb  -  0 ,028$h,  z  -  0 . 7h) 


i ho  phase  difference  of  each  elementary  wave 
does  not  change  bow  wave  resistance  from  the 
linear  theory. 

The  stern  wave  amplitude  is  exactly  the 
same  as  the  Linear  value 

A  (“)  -  P  +  LQ 

s  s  .s 

hut  the  phase  d  i  t  to rune l  (••)  should  bo  corn- 

put  oil  by  tho  ray  theory  together  with  the  ray 
path  which  U  shown  in  Figure  U.  Then  It  is 
also  obvious  that  the  stern  wave  resistance  is 
the  same  as  the  linear  stern  wave  resistance. 
The  total  wave  resistance  may  be  obtained  simi- 
lai ly  by  considering  that  the  total  wave  ampl i- 
t  udo  is 


considered  to  be  totally  missing  in  too  wave 
resistance,  l  in  the  Srettensky  formula  for  how 
waves  would  start  not  from  zero  but  from  a 
certain  number  minimum  .1  *  dj  such  that 


whore  "u.j  tho  value  of  -\t  tram  which  bow 
waves  start  to  reflect.  If  ow,  *  u,,  t.  -  1, 
and  £.  •  2  foj  .1  •  J  i  ,  and  If  o  * 

«  j  -  Z  for  ,1  >  Jp 

The  result  is  shown  in  Figure  12  where  a 
considerable  shift  of  the  phase  of  hump  and 
hollow  of  the  wave  resistance  due  to  bow 
stern  wave  interaction  is  noticeable. 
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Hoii-,  the  bow  and  stern  wave  mti'uu  t  ion  ap¬ 
pears  in  the  wave  resistance.  That  is,  onlv 
the  int  or  ai  t  ion  term  changes  due  to  tin*  phase 
change  caused  by  the  nonuniform  flow.  Ibis 
luct  » ex.it  i  1  v  the  same  as  in  two-dimensional 
t  heui  y  . 


The  actual  computation  of  wave  resistance 
is  performed  by  the  Srettensky  formula  using 
the  relation 


-  cos'1  (— ->  -  un  '  id  ‘  -  I)'5  (45) 

and  the  corresponding  values  ot  sj  are  obtained 
hv  interpolat  Ion .  When  a  portion  of  elementary 
waves  neat  a  0  is  reflected  1  tom  the  ship 
hull .  the  large*  part  of  the  energy  in  this 
portion  of  element arv  waves  will  be  dissipated 
hv  breaking  waves,  and  the  wave  resistance  till 
decrease,  but  the  momentum  loss  due  to  breaking 
waves  will  inctease.  l!  such  energy  was  con- 
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Figure  11  -  Kay  Paths  of  Stern  Waves  of 
a  Wigley  Hull  (b  -  0.2,  h  »  0.0625) 


Figure  12  -  Wave  Resistance  of  the 
Wigiey  Ship  (b  »  0.2,  h  *  0.0625) 

Discuss  ion 

Recently  Kggcrs^  obtained  a  dispersion 
l  elation  tiom  a  low  speed  t  rev-sur  r  ac  e  boun¬ 
dary  condition  which  was  .«  slight lv  modified 
version  ol  the  one  huba***  used.  Rogers  showed 
that  there  was  a  sn-.l  l  region  near  the  stag¬ 
nation  point  where  the  wave  number  became 
negative;  since  this  was  not  permissible  lor  a 
wave  and  it  could  be  mterpteted  to  mean  that 
there  was  no  wave  in  the  region.  He  suggested 
that  his  torm  might  help  alleviate  the  sen- 
iltivitv  ol  l lie  initial  condition  to  the  ray 
paths.  iV.  voiding  to  the  energy  conservation 
law**  ol  the  wave  propagating  through  non- 
uniform  Hew  the  wave  energy  Mux  is  proportion¬ 
al  to  the  wave  numiict  along  tin  ray  tube. 

When  the  wave  number  along  a  ray  is  considered » 
although  the  Keller  dispersion  relation  has  an 
iniimte  wave  number  onlv  at  the  stagnation 
point  ,  the  I  >;ge  t  :•  dispersion  relation  has  an 
mi  mite  wave  number  in  •  he  flow  near  the  bow. 
I'lins,  neat  such  a  singular  paint  or  a  singular 
line,  waves  may  break  and  the  present  theory 
cannot  be  applied.  When  tggers*  equation  was 
nu orporated  into  the  rav  equation  in  the 
ptesont  lay  computer  urogram  it  was  t ound  that 
the  tav  path  was  still  verv  sensitive  to  a 
small  change  et  initial  value  lx,y,  )  and  the 
.  urvo-l  tav  *»iiM  reflected  iron  the  ship  hull. 


In  the  experiments'  conducted  at  Tokyo 
l  tuvers  U  v ,  the  second  caustic  can  be  noticed, 
and  near  the  second  caustic  the  f low  field  is 
violently  different  from  the  linear  theory. 
However,  it  seems  that  extreme  caution  is 
needed  to  distinguish  the  second  caustic  from 
the  first  caustic.  The  case  of  experimentsJ 
with  wedges  is  interesting  because,  according 
to  the  present  theory,  there  cannot  be  a  second 
v.iustu  on  the  wedge  although  the  tiiSv  caustic 
should  be  there.  However,  the  wave  phases  oi 
the  wedge  should  be  advanced  and  quite  sensitive 
to  the  draft*  and  beam- length  ratios  due  to  the 
nonuni form  flow  caused  by  the  wedge.  On  the 
other  hand,  the  present  theory  is  still  not 
exact  although  if  takes  into  account  the  effect 
on  tin*  propagat ion  of  water  waves  ol  nonunitorm 


flow  due  to  th«  double  model.  Thus,  the  non¬ 
linear  effect  of  water  waves  is  not  totally 
analyzed  here.  Nevertheless,  the  present 
theory  supplies  a  great  deal  of  hope  to  an 
entiiely  different  ap  roach  to  the  ship  wave 
theory  -  the  ray  theory. 

Although  the  strictly  linear  amplitude 
L unction  is  used  for  simplicity  In  the  present 
study,  a  slightly  improved  amplitude  function 
may  be  easily  incorporated  by  adding  the 
sheltering  effect  or  other  effects.  However, 
it  should  be  noted  here  that,  In  any  case,  the 
bow  wave  resistance  and  the  stern  wave  resis¬ 
tance  are  th*  same  as  the  results  without  the 
ray  theory  and  the  effect  of  the  ray  theory 
would  appear  as  a  shift  of  the  phase  of  hump 
and  hollow  in  the  total  wave  resistance.  There¬ 
fore,  the  computational  results  would  not  match 
the  results  of  towing  tank  experiments  unless 
the  viscous  boundary  layer  effect  on  the^tern 
wave  or  some  other  effect  is  considered. 
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Appendix  A 


For  the  computation  of  ray  paths  of  a 
ship,  the  flow  velocity  and  its  derivatives  on 
(x.y.o),  u ,  v  ,ux ,  v  »Uy  are  needed.  A  Wig’.ey 
hull  has  the  double  model  source  distribution 


m  *  b  ( 


z.  2. 

-2xt  f  1)  1  -  (jji.)  I 


0  <  x  <  I.  y  «  0,  h  >  t.  >  -h 


Thus , 


h  l 


-u  (x,y,o) 


2  /  /  m  &  $  d*i  dD  ' 


0  o 
h  l 


I  /  m~M 71  d*'  dz'  *  1 


'I  l 


».  1  -  — 


2b 


i  [ _ u 

'  J  r  <*!  ■ 


iTd*i 


h  j  -  — L_ 

1  2 

h“  dz 


It 


(x |  *  o) 


1  -  1 


+  4b  £  b  log  j  Xj  - 
+  r(Zj  »  h)}  - 


13 


Uj  -  x)  y 


— — rr~  d,i 

(7^  +  y  )  r 

J  U1  '  X)  d«, 


I  4  (x  -  x) 

y  _1 - - 

(*i  +  y2>  r 


‘l  -  X)  (Z1  ~  y2)>dz.  1  1 

r  Jx,-0 


c.y.O)  -  2  f  /»% 


dXj  dzj 


(  2x2  - 

-2  by  — - 2 

i  L<*  +  y 


2x  -)xtl _ 

+  y“)  ru^  “  l) 


(l  ♦  V  -  — — -^7 

h  r  (Xj  “  1)  !• 


r(Xj*l)  h 


2  +  r(Xj-l) 


2x  -x  +  2z 


227  .! 

-J- - dz! 

h‘  r(x,  *0) J 


r(x^*0)  h 


1  r  (x  |  - 1 , 

r,-0) 

-u 

X 

-4b  1015  |  r  (X|*n, 

i:i"0) 

4 

i-  (>j  r(x,’l) 

_  *_± 
2 

log  i  z,  -  rU,-l> 

(x-1) 

')  2 

al 

+  2  | 

+  r  1  i  j 

1l 

rU.-l)| 

y*  *i  r(v°> 


log  lij  +  r(Xj-O)  I 


2  {(  —  +—)  l°g  + 

(V  2  4 

1  h 

r(x  -0)  | - -  r(x.”0)j 

1  *  'J, 


2  v?  2 

where  •»  (x-l>  ’  +  y 


2 

-u  *  1  j  j  "  IJsJ  (7)  dxi  d*i 

o  o 

2by  [,,,  +  y2}  td.W 


l  .  l 
+  r(l  ,h)  h  "*  r(0,h)  h 


1  ,  (h  +  r(l,h){h  +  r(o.h)) 

- ;log  — — - - — 

h2  r(l,0)  *  r(O.O) 

_ J. - 1 

(x2+y2/  r(o,h)J 

2 

r  (1-x)  (l  +  Mdzl  + 

+4by  {  I _ v  h _ 

0  (*i  +  y2)  r(Vl> 

}  "('*?)  „ 

+  7~7  TT  dzl 

l  (Z1  ♦  v‘)r  (xr0) 


-i,  /  JZ- 

h‘  2  rfXj-l) 


-  .2  )  ri^o)  diil 


v 


b(by2  +  )* 


X 


4 


) 

t'tj-x)*  - 
-  o 

iX|-x)  y‘ 


>f 


o 


zj>y  j  r 


2.  culler,  J . ,  "The  Ray  Theory  of  Ship  Waves 
and  the  Class  of  Streamlined  Ships/'  J.  of 
Fluid  Mechanics,  Vol.  91,  Part  3,  pp,  465- 
488  (1978). 

3.  Miyata,  H.,  "Characteristics  oi  Nonlinear 
Waves  in  the  Near-Field  of  Ships  and  Their 
Effects  on  Resistance,"  Proe.  of  l 3th  ONK 
Symposium  on  Naval  Hydrodynamics  (1980). 

4.  Uracil |  F.,  "Steady  Wave  Patterns  on  a  Non- 
Uniform  Steady  Fluid  Flow,"  J.  cf  Fluid 
Mechanics,  Vol  9,  pp.  333-346  (1960). 


2  _ h _ 

3  [(Xj-x)‘  +  y"  r  U  j  *h) 


Shen  M.C.,  K.E.  Meyer  at:d  J.B,  Keller, 
"Spectra  of  Water  Waves  in  Channels  and 
Around  Islands,"  The  Fhyslcs  of  Fluids, 
Vol  11.  No.  11,  pp.  2289-2304  (1968), 


+  -jU-2x) 


-h 


(x 


2  21 

-x)  (x-Xj)  +y  JrlZj-h) 


11  +  <zi~h> 


(h“  +  y")  (Xj-x)  y 


Xj-0 


,  j  2y"  +  !  -2x)  (x  -x)  h 

6  b>-‘  ' - 2 - i - - 

ih  \  (Xj-x)  4  y  rUj*h) 


6.  Yim,  B.,  "Notes  Concerning  the  Ray  Theory 
of  Ship  Waves,"  Proc .  of  the  Continued 
Workshop  on  Ship  Wave-Res iatance  Comput- 
tations,  pp,  139-153  (1980). 

7.  Inui,  T.,  H.  Kajitani,  H.  Miyaca,  M, 
Tsuruoha,  and  A.  Suzuki,  "Nonlinetr  Proper¬ 
ties  of  Wave  Hiking  Resistance  of  Wide -Beam 
Ships,"  J .  of  the  Society  of  Naval  Arch,  of 
Japan,  Vol.  146,  pp.  )8-2&.»(Dec  1979). 

8.  Lighthill,  J.,  "Waves  in  Fluids,"  Cambridge 
University  Press,  Cambridge,  London  pp.  332- 
337  (1978). 


+  2  lo 


h  +rUj“h)  2,i  \ 

L>r  Uj-0)  rCaj-hy 

h  . 

J  -At 

J  (Z7  +  V 
O  I 


2  (x  j - x)  ( l-2x) 
3hJ 


)  r 


x(‘:0 


9.  Doctors,  L.  and  0,  Dugan,  "Comparison  of 
Nonlinear  Wave- Resistance  Theories  for  a 
Two-Dimension «l  Pressure  Distribution," 

.1.  of  Fluid  Mechanics,  Vol.  98,  Part.  3, 
pp.  647-672  (1980). 

10.  ^im,  B.,  "Design  of  Bulbous  Bow  by  the  Rdy 
Theory  of  Ship  Waves,"  to  be  published. 


v  (x , y ,0) 

y 

vneit  iui»b)  r(Xj»a,  4j=b) 

.n  these  'xpressions  the  integrals 


11.  Yim,  B.,  "A  Simple  Design  Theory  and  Method 
for  Bulbous  Bows  of  Ships,"  J.  of  Ship 
Research,  Vol.  18,  No,  3,  pp.  141-152 

(Sep  1974). 

12.  Havelock,  T.H.,  "Wave  Patterns  and  Wave 
Resistance ,"  Tran,  of  the  Inst,  of  Naval 
Architects,  Vol.  76,  pp.  430-443  (19J4), 


n 


J  2  2 

4  y  •  r 

20 

•ii'i  available  in  closed  forms. 

References 

Inui,  T.  and  H.  Kaiitam,  "Study  on  Local 
Non-Linear  Free  Surface  Effects  in  Ship 
Waves  and  Wave  Resistance,"  Schif f stechnik , 
Band  24,  Heft  118  (Nov  1977). 


13.  Inui,  T.,  "Wave -Making  Resistance  of 
Ships,"  Trans,  of  SNAME,  Vol.  70,  pp.  283- 
31  1  (1962) . 

14.  Lamb,  H.,  "Hydrodynamics,"  Dover  Pub., 

New  York.  N.Y.,  p.  395,  (1945). 

15.  Havelock,  T.H.,  "The  Calculation  of  Wave 
Resistance,1'  Proc,  of  the  Royal  Society 
(London),  Vol.  A144,  pp.  514-521  (1934). 

16.  Srettensky,  L.N.,  "On  The  Wave  Making 
Resistance  of  a  Ship  Moving  Along  in  a 
Canal,"  Ph i losophica 1  Magazine,  Vol.  22, 

7th  Series  (1936). 

17.  Eggers,  K. ,  "On  the  Dispersion  Relation 

and  Exponential  Variation  of  Wave  Components 


15 


satisfying  the  Slow  Ship  Oil' feron'  ini 
Equation  on  the  I'ndlslurbed  Free  Surface," 
International  .hunt  Research  on  S»  udv  on 
Uh.i1  Nonlinear  Effect  in  Ship  Waves, 
Research  Report  1979,  Edited  by  T,  Inui, 
pp,  4  )-bJ  l Apr  1980). 

Hob. i ,  I.,  and  M.  Hava,  "Numerical  Evaluation 
of  Wav.*  Resistance  Theory  lor  Slow  Ships," 
Pro-- ,  Jnd  International  Conference  on 
Numerical  Ship  Hydrodynamics,  ITuversiiv 
of  California,  Berkeley,  pp.  1  7-J9  (1977)  . 

Moteno,  M.,  1.,  Perez-Ro  pis  and  !,  l.andweber 
"Effect  of  Wake  on  Wave  Resistance  of  a 
Ship  Monel,"  Iowa  Institute  of  Hydraulic 
Reseairh,  Report  180  (Aug  197*>)  . 


Yim,  B.,  "Theory  of  Ventilating  on 
Cavuating  Flows  about  Symmetric  Surface 
Piercing  Struts,"  David  W.  Taylor  Naval 
Ship  Research  and  Development  Center, 
Report  4b  lb  l Sep  l 9 7b). 


NITIAL  DISTRIBUTION 


Copies 

1  CHONR/438 

2  NRL 

1  Code  2027 
1  Code  2627 

4  USNA 

1  Tech  Lib 
1  Nav  Sys  Eng  Dept 
1  Bhattachcryya 
1  Calisal 

2  NAVPGSCOL 

1  Library 
1  Garrison 

1  NELC/Lib 

2  NOSC 

1  Library 
1  Higdon 

1  NCEl./Code  131 

13  NAVSEA 

1  SEA  031,  R.  Johnson 

1  SEA  031,  G.  Kerr 

1  SEA  031,  C.  Kennel 

1  SEA  03R,  L.  Benen 

1  SEA  03R,  Dilts 
1  SEA  03R,  N.  Kobitz 

1  SEA  03R,  J.  Schuler 

1  SEA  312,  P. A.  Gale 

1  SEA  312,  J.W.  Kehoe 

1  SEA  321,  E.N.  Comstock 

1  SEA  321,  R.G.  Keane,  Jr. 

1  SEA  61433,  F.  Prout 

1  PMS  383,  Chatterton 

1  NADC 

12  DTIC 

1  NSF/Engineer ing  Lib 

1  DOT/Lib  TAD-491.1 

1  NBS/Klebanoff 


Copies 

1  MARAD/Lib 

4  U.  of  Cal/Dept  Naval  Arch, 

Berkeley 
1  Eng  Library 
1  Webster 
1  Paulling 
1  Wehausen 

2  U.  of  Cal,  San  Diego 

1  A.T.  Ellis 
1  Scripps  Inst  ib 

2  CIT 

1  Aero  Lib 
1  T.Y.  Wu 

1  Catholic  U.  of  Amer/Civil  &  Meeh 

Eng 

1  Colorado  State  U./Eng  Res  Cen 

1  Florida  Atlantic.  U. 

1  Tech  Lib 

1  George  Washington  U./Feir 

1  U.  of  Hawaii/St.  Denis 

1  Harvard  U. /McKay  Lib 

1  U.  of  Illinois/ J .  Robertson 

2  U.  of  Io>»a 

1  Library 
1  Landweber 

1  Johns  Hopkins  U. /Phillips 

1  U.  of  Kansas/Civil  Eng  Lib 

1  Lehigh  U. /Fritz  Eng  Lab  Lib 

1  U.  of  Maryland/Dept  of  Aerospace 

Eng,  Plotkin 


17 


SNAME/Teeh  Lib 


Copies 

5 


3 


1 


A 


1 

2 


2 


3 


1 


MIT 

1  Young 
1  Mandel 
1  Abkowitz 
1  Newman 
1  Sclavounos 

U.  of  Mich/NAME 
1  Library 
1  Beck 
1  Ogilvie 


Copies 
1 

1  Bethlehem  Steel/Sparrows  Point 
■  1  Bethlehem  Steel/New  York/Lib 

1  Exxon,  NY/Design  Div,  Tank  Dept 

1  Genet al  Dynamics,  EB/Boatwr ight 

1  Gibbs  &  Cox/Tech  Info 


U.  of  Notre  Dame 
1  Eng  Lib 

New  York  U./Courant  Inst 
1  A.  Peters 
1  J.  Stoker 

SIT 

1  Breslin 
1  Savitsky 
1  Dalzell 
1  Kim 


A  Hydro  nau  t  ic.s 

1  Library 
1  Elsenberg 
1  Tulin 
l  Hsu 

1  Lockheed,  Sunnyvale 

1  Potash 

1  Newport  News  Shiphu  1  Iding/I.ib 

1  Oceanic;; 


U.  of  Texas/Arl  Lib  1 

Southwest  Rest  Inst  1 

1  Applied  Mech  Rev 
1  Abramson 

Stanford  Res  Inst 
1  Civ  Engr  Lib 
1  Street 

VPI 

1  Nayfeh 
1  Mook 

U.  of  Washington 
1  Eng  Lib 
1  Mech  Eng/Adee 

Webb  Inst 
1  Library 
1  Lewis 
1  Ward 

Woods  Hole/Ocean  Eng 


Sperry  Rand/Toch  Lib 

Sun  Shipbuilding/Chief  Naval  Arch 

American  Bureau  of  Shipping 
1  Lib 
1  Cheng 

Maritime  Research  Information 
Serv ice 

CENTER  DISTRIBUTION 


Code 

Name 

0120 

D,  Jewel  1 

1  102 

C.D.  Elmer 

11/ 

R.M.  Stevens 

1170 

1170 

O.R.  Lamb 

S.  Hawkins 

1  300 

W.R.  Morgan 

2 

1 

Copies 

1 

l 

1 

1 

1 


18 


CENTER  DISTRIBUTION  (Continued) 


Copies 

Code 

Name 

1 

1504 

V..T.  Monecelle 

1 

1506 

D.  Ciealowski 

1 

1520 

W.C.  Lin 

1 

1521 

W.G.  Day 

1 

1521 

A.M.  Reed 

1 

1522 

G.F.  Dobay 

1 

1540 

J.H.  McCarthy 

10 

1540 

B.  Yirn 

1 

1540 

T.T.  Huang 

1 

1542 

C.M.  Lee 

1 

1542 

M.  Chang 

1 

1552 

M.B.  Wilson 

1 

1560 

Division  Head 

1 

1561 

G . C .  Cox 

1 

1561 

S.L.  Bales 

1 

1562 

D.D.  Moran 

1 

1562 

E.E.  Zarnick 

1 

1562 

Y.S.  Hong 

1 

1563 

W.E.  Smith 

1 

1564 

J.P.  Feldman 

30 

5211.1 

Reports  Distribution 

I 

522.1 

Unclassified  Lib  (C) 

522.2  Unclassified  Lib  (A) 


19 


